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^ ■ Abstract 



We study the nuclear (A) dependence of the European Muon Collaboration 
(EMC) effect at high values of x (x > 0.6). Our approach makes use of con- 
ventional nuclear degrees of freedom within the Relativistic Impulse Approx- 
imation. By performing a non-relativistic series expansion we demonstrate 
that relativistic corrections make a substantial contribution to the effect at 
x > 0.6 and show that the ratio of neutron to proton structure functions ex- 
tracted from a global fit to all nuclei is not inconsistent with values obtained 
from the deuteron. 
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The famous high energy deep inelastic lepton scattering results obtained in 1983 at CERN 
by the European Muon Collaboration (EMC) jl]] and at SLAC in 1984 (now known as the 
EMC effect) showed that quark momentum distributions are modified in the nuclear medium. 
The scattering is well described by assuming that the scattered lepton interacts with a bound 
quark by exchanging a virtual quanta with four-momentum Q 2 = q 2 — u 2 , where q = |q| and 
v are the three-momentum and energy carried by the quanta, respectively. The inclusive 
cross section depends on the Bjorken variable x = Q 2 /(2Mnu), which is identified with the 
fraction of the total longitudinal momentum carried by the struck quark. In this Letter 
we will discuss recent progress in the interpretation of the data for large quark momentum 
fractions x > 0.6. 

Figure 1 shows the ratio Ra(x) of the nuclear structure function per nucleon, F^(x)/A, 
to the deuteron structure function per nucleon, F 2 D (x)/2. The data are for an iron nucleus, 
but the results are similar for all nuclei with mass number A > 3. If the nuclear medium 
had no effect on the quark momentum distribution, the ratio would be unity; the (up to 
20%) deviations of the nuclear structure functions are direct evidence for the effect of the 
nuclear medium. 

Throughout the years, a variety of models have been proposed to explain the EMC effect 
at large x (for a review see |f|). Some invoke additional (sometimes exotic) components of 
the nuclear wave function, while other, more conventional models describe the EMC effect in 
terms of nucleon binding. In this latter picture, the rise above one as x — > 1 is due to Fermi 
motion. In Ref. Q, referred to as CL in this paper, it was pointed out that a nonrelativistic 
calculation based on binding plus short range nucleon-nucleon (NN) correlations (generated 
by the NN potential) were sufficient to account for most of the effect at x < 0.5, but the 
result for large x (the long dashed line in Fig. 1) was much too small. As this model requires 
only a very small number of parameters which can be determined to high precision from 
other data ||, agreement cannot be obtained by adjusting the parameters, leaving open the 
possibility that exotic components could play a role in the explanation of the effect ||. 

In Ref. [0, which we will refer to as GL, we introduced a new relativistic formulation of 
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the impulse approximation (which we refer to as the RIA), based on the relativistic spectator 
model || . Using this formulation, we found that we can explain the EMC effect in the region 
of large x entirely in terms of conventional nuclear degrees of freedom, as shown in Fig. 1. 

In this Letter we reinforce these conclusions by showing that (i) the theoretical uncer- 
tainty of the RIA calculation is not large, so that the agreement presented in GL is no 
accident, (ii) the relativistic effects can be rather simply parameterized, and hence their 
physical origin easily understood, and (iii) the new relativistic theory removes a descrep- 
ancy between the results obtained for the neutron structure function, F2 n (x), in two different 
ways: from measurements of the deuteron, and from global fits to nuclei with all values of 
A. 

To begin the discussion, we recall that the A dependence of the experimental ratio R A (x) 
has been generally parametrized as a product of a function of x times a function of A: 

R A (x) - 1 ~ a(x)P(A) . (1) 

This is referred to as "factorization". Note that the parameterization || Ra{%) oc A a ^ takes 
this form if the exponent a(x) is sufficiently small. A factorized form for R A {x) can also 
be obtained from a naive nonrelativistic IA description. In this case the structure function, 
F 2 A , is given by a linear convolution formula PlJlOll 

F 2 A (x) 

where F 2 is the structure function for an off-shell nucleon (which is assumed to have the 
same form as the on-shell one), z = A(k ■ q)/(PA ■ q) — k + /M^ is the light-cone momentum 
fraction of the struck nucleon with four-momentum k = (k , k), and Ma are the masses 
of the nucleon and nucleus A, respectively, and Ja{ z ) is the nucleon light cone momentum 
distribution. 

To obtain the factorized form ([I]) from the convolution formula (0), it is useful to exploit 
the fact that the nuclear momentum distribution /a(^) is sharply peaked around z — 1, and 
expand the factor F 2 (x/z) in Eq. (§) in powers of (1 — z) around z = 1 
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The coefficients of the expansion, accurate to order 1/Mjv, are therefore 
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where (T) is the average kinetic energy of the nucleon in the nucleus and (E) = (Ma-i) + 
M/v — M A is the average removal energy, with (M A -x) the average mass of the spectator 
A — 1 nuclear system (in this discussion we neglect the recoil energy of the A — 1 system). 
Details of the derivation of these coefficients are discussed in CL. The coefficient c is just 
the normalization of the light cone momentum distribution, and c% and c<i can be related to 
(T) and (E) by exploiting the connection between f A (z) and the nucleon three- momentum 
distribution, n A {k) 

(E) k {{ \ 
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where k — |k| is the magnitude of the three momentum of the struck nucleon, and k,, its 
component in the direction of the q. 

Note the presence of the the factor z (sometimes referred to as the flux factor) in these 
equations. This quantity was omitted from some early papers on nuclear deep inelastic 
scattering because incorrect assumptions were made in connecting the relativistic formalism 
with the nonrelativistic distributions actually used in the calculations. Its effect on nuclear 
structure functions was emphasized in [O, but here we wish to emphasize that the flux factor 
does not change the normalization by more than a few percent, which is not numerically 
significant in the discussion of the EMC effect at large x, and it has no effect on the coefficient 
C2- The principal effect of the flux factor is to add the term — 2(T) / (3Mjy) to the c\ coefficient 



in Eq. (|J). This has a pronounced effect, reducing the size of this coefficient by almost a 
factor of two, and decreasing the size of the EMC effect predicted by the nonrelativistic 
impulse approximation (for a detailed discussion, see CL). 

A factorized equation of the form (|1|) can be obtained by using the energy-weighted sum 
rule [|IJ 

(E) = ^(T)+2e^ 3 -(T), (6) 

where e is the binding energy per nucleon for nucleus A, and the second expression uses the 
approximation 2e « (T) /2. With these approximations we obtain 

«(*) = r^j^ + r 2 ^f > m = f)/m n . (7) 

Therefore, in nonrelativistic IA, the A and x dependencies of Ra — 1 factor, with the A 
dependence given entirely by the average kinetic energy of a nucleon in a nucleus, and 
the x dependence contained entirely in the term which depends on derivatives of the free 
structure function F 2 . In CL it was shown that realistic momentum distributions including 
NN correlations yield values of (T) which are large enough to reproduce the EMC effect at 
x < 0.5, but the expansion (0) does not explain the behavior of Ra at higher x (as already 
shown in Fig. 1). 

Calculations performed with the Relativistic Impulse Approximation (RIA) introduce 
some important corrections to the IA convoution formula, Eq. @. In RIA, the nuclear 
structure function becomes 



F 2 A (x) 



RIA 



dzf* IA (z,x), (8) 



with 
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f% IA (z,x) = 2nM N z / dkkn A (k)F 2 N (y,k), (9) 

where y = rjx/z (rj = AMn/Ma — 1) is the momentum fraction carried by a quark inside 
the bound nucleon, F 2 N (y, k) is the structure function for the bound (off-shell) nucleon which 



depends on the longitudinal and transverse momentum of the nucleon through z and k, and 
riA{k) can be related to a covariant nuclear spectral function (which is not known, but is 
determined by its relation to riA{k)). Further analysis of the relativistic kinematics permits 
us to express the off-shell nucleon structure function, F 2 N (y, k), as a product of a relativistic 
phase space factor, P(y,y') times the on-shell structure function of a shifted argument, 
F 2 N (y'), 

F 2 N (y,k)=P(y,y')F 2 N (y') (10) 

where P(y, y') describs the phase space of the spectator quarks (and satisfies the condition 
P{ViV) = 1)) an d y' is the value of y shifted by the relativistic kinematics: 

y' = y(l- |A) + | ^(b 2 (y) + lyA) 2 + y{\ - y)A - %)} , (11a) 

1 ml/Ml 1 , . , . 

b(y) = -^-^--(l-y) (lib) 

A _ m 2 -kl _^ (M A - M A -if n M A (Ea-i ~ M A -i) ( . 

~ M% ~ Ml +Z Ml ' (UCj 

where mx > Mjy is the mass of the spectator quarks with relativistic phase space P(y,y'), 
m is the mass of the struck quark, k^ = k^ — k 2 ^ m 2 is the square of the four-momentum of 
the off-mass-shell struck quark, and Ma-i is the mass of the recoiling A — 1-nuclear system. 
Note that A is a measure of how far the struck quark is off-shell, and is zero if it is on-shell. 

This important connection between the off-shell and on-shell nucleon structure function 
was first obtained in GL. The most significant difference between the non-relativistic convo- 
lution formula, Eq. @, and its relativistic counterpart, Eq. (||), is the explicit dependence of 
the nucleon structure function, F 2 N , on A; when A = one can easily see that y' = y and the 
IA convolution formula, Eq. @, is recovered. We emphasize that Eqs. (I9Q— (|llc|) have been 
derived from considerations of relativistic kinematics only; any dynamical dependence of F 2 
on the nuclear medium, i.e. an explicit dependence on k 2 not due to relativistic kinematics, 
has been disregarded. This means that our relativistic formulae are truly consequences of 
the assumption that the nucleons are not modified by the nuclear medium. 
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Both the A and x dependence of the EMC-effect at large x is significantly altered by 
the parameter A, which includes the relativistic corrections to the IA. We can easily display 
(approximately) the effect of these relativistic corrections by expanding f A IA (z, x) in powers 
of the small quantity A. We obtain: 

ff A {z,x) ~ f A (z)F 2 N (y) + (A) ± G N (y) + 0(A 2 ) , (12) 

where 



A=0 



dy' 



° F 2 N (y')P(y,y>) 



(13) 
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and 
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(A) ± = 2nM N z dkkn A (k) A (14) 

is the average value of A over the nucleon transverse momentum k. Note that the first term 
in Eq. (|i~2D (the one proportional to f A ) is identical to the nonrelativistic IA result given in 
Eq.©. 

Expanding the first term in Eq. fll~2|) around z = 1, as we did before, gives 



R A (x) - 1 ~ ai {x)^{A) + a 2 (x)(3 2 (A) , (15) 



where ai(x) and fl\(A) are the a(x) and (3(A) given in Eq. (0), and 
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p 2 (A)= (A(z)) = 2irM N dz z dkkn A (k)A. (16) 
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The new, relativistic correction term /?2(^4) is the value of A averaged over both longitudinal 
and transverse momentum variables. 

To summarize: the RIA still gives a result in which the x dependence of the nuclear 
structure function is rescaled by the motion of the nucleons, and the A-dependence is still 
governed by average properties of nucleon dynamics which can be readily calculated with 
high accuracy by present day nuclear models. The important difference in the region of 



large x is that the A and x dependencies of the EMC effect cannot be written as a product 
of a single factor of a function of x times a function of A, but require the sum of two such 
products. This is a consequence of relativistic effects, which give the second term in Eq. (fTBD . 

We now turn to a discussion of the three points mentioned at the beginning of this 
letter. The shaded area in Fig. 1 are the predictions of the RIA, including an estimate of the 
theoretical error. This error includes uncertainties due to different parametrizations of the 
free nuclear structure functions (shown by the three dotted lines in the figure) and variations 
in the nuclear parameters of our theory, namely the average removal energy (E), and the 
value of the mass parameter mx for the spectator quarks. Of these, the only significant 
error comes from the dependence of the theory on the (unknown) mass of the spectator 
quarks (which is expected to be a mass close to, but larger than, M N ). The two solid lines 
shown in Fig. 1 correspond to nix — 940 and 1800 MeV, showing that the predictions are 
insensitive to the preceise value of this parameter. We conclude that the agreement between 
theory and experiment is no accident. 

Finally, we turn to the question of the extraction of the neutron structure function from 
experimental data. The usual way to obtain the neutron structure function is to measure 
the ratio of the deuteron to proton structure functions. Ignoring nuclear motion effects in 
the deuteron, this ratio is 

^l-l = ^l = R(x) (17) 

However, if the nucleon structure functions are not modified by the nuclear medium, and 
if our theory of the EMC effect is good enough, then the ratio R(x) could be, in principle, 
also extracted from measurements on any other nucleus, even if the nuclear recoil effects 
are large for that nucleus. The value of R(x) obtained in this way should agree with the 
result obtained from the deuteron. This idea is best tested by extracting R(x) from a global 
fit to all available nuclei (A = 2,4,9,12,27,40,56,108,197) |12|,|13| using the theory only 
to specify the A dependence through the functions /3%(A) and ^(A). One can then not 
only extract R(x) from the global fit, but also the "coefficients" a\{x) and a2(x), and hence 
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determine all of the x dependence of the EMC effect in an independent way which is not 



biased by the theory This idea has already been carried out in Refs. [[L2|,|n|, where the 
A-dependence was modeled according to an empirical formula [TO] . 



The results of this analysis ||12| , P~3f gave rise to a puzzling situation: while at 0.3 < x < 0.6 
the extracted values of R were in good agreement with the R extracted from deuteron, at 
higher x (x > 0.7) there appeared to be a significant discrepancy. In particular, in the limit 
x — > 1, the deuteron data seems to extrapolate to the lower bound of 1/4, while the global 
fit extrapolates to a much higher value (compatible with the 3/7 suggested by Farrar and 
Jackson @). 

We have extracted the value of R from a global fit to the EMC data, using the A- 
dependence suggested by Eq. (|1^). The ratio is obtained by fitting the formula 

R *ex P= (T* A exp {x) _ Na'Kx) + Za v A {x) _ 



N 

-R(x) + 



N 

1 + -R[x) 



[a 1 {x)i3 1 {A) + a 2 {x)(3 2 {A)} , (18) 



where the first line defines the experimental ratio R*^ xp []13|] , and (3i(A) and fl 2 (A) are the 
theoretical input. The fit determines R{x), aci(x), and a 2 (x). The results of the fit, shown 
as black squares in Fig. 2, are compared with the theoretical results (solid curves) obtained 
from the parameterizations of F 2p and F 2n (obtained from deuteron data). 

The open squares in Fig. la were obtained in Ref. ]ll| (referred to as BDR). Note that 
the trend of the BDR fit clearly suggests disagreement with the deuteron data (the solid 
line) and also suggests that R — > 3/7 as x — >1, while our fit is full agreement with the 
deuteon data, and extrapolates to a limit near 1/4 as x — > 1. 

It is now becoming increasingly clear that the nature of the EMC-effect cannot be unrav- 
elled unambigously without addressing the problems inherent in the formalism (relativistic 
covariance, meson degrees of freedom and final state interactions) [[HJ. Our formula includes 
a breakdown of factorization due to relativistic effects. A margin of accuracy of the theory 
must be anyway specified and compared to the magnitude of the effect. We therefore present 
our results along with the theoretical "error". More accurate extractions could be obtained 
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by measuring the EMC-effect at large x for a wider number of nuclei. 
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FIGURES 

FIG. 1. The EMC data for 56 Fe. The dashed line is the nonrelativistic calculation of CL. The 
shaded area is the relativistic calculation of GL, including an estimate of theoretical uncertainties. 

FIG. 2. Results of a global fit to the EMC data. The solid squares (with error bars) are 
obtained from a three parameter fit to Eq. (|l8|); the open squares (displaced slightly in x so that 
they can be separated from the solid squares) are a two parameter fit with the relativistic term 02 
omitted. The short bars in (a) are the limits 3/7 and 1/4. The solid lines are discussed in the text. 
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